Quantum analogue of the spin-flop transition for a spin pair 
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Quantum (step-like) magnetization curves are studies for a spin pair with antiferromagnetic cou- 
pling in the presence of a magnetic field parallel to the easy axis of the magnetic anisotropy. The 
consideration is done both analytically and numerically for a wide range of the anisotropy constants 
and spins up to S > 100. Depending on the origin of the anisotropy (exchange or single-ion), the 
magnetization curve can demonstrate the jumps more than unity and the concentration of the unit 
jumps in a narrow range of the field. We also point the region of the problem parameters, where the 
behavior is quasiclassical for S = 5, and where system is substantially quantum in the limit S — » oo. 
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Mesoscopic magnetic systems whose total spin S is 
large, but that exhibit quantum effects associated with a 
finite value of the spin, are have been studied both experi- 
mentally an theoretically in the last decade. For instance, 
such systems are high-spin molecules with S ~ 10, mag- 
netic clusters with S ~ 100, and magnetic dots (small 
magnetic particles with size 50 - 100 nm), see Refi£. 
Objects that demonstrate quantum properties in macro- 
scopic (more accurate, mesoscopic) scales are important 
for the physics of the magnetism and its applications: the 
implementation of the idea of a quantum computes^ and 
information recording devices i The simplest quantum ef- 
fect that is manifested in such systems is associated with 
the quantization of the projection of the total spin. It 
consists in a step-like change of the magnetization when 
the external magnetic field changes continuouslyiS Such 
a behavior is observed for many systems with an anti- 
ferromagnetic (AFM) interaction, see for a review^ The 
most studied systems are spin dimers with AFM interac- 
tion, for instance, pairs of high-spin molecules Mn4 with 
the maximum total spin 2S* = 9, seeSii Measurements 
are also carried out for spin triplets, quartets, etc££ The 
analysis of jumps on the magnetization curves is useful 
for experimental determination of material constants. 5 

For realization of this method, it is important to 
note that if anisotropy terms are absent in the Hamil- 
tonian QJ, i.e. B = and k = 0, the problem has an 
exact solution for an arbitrary value of the spin. The 
eigenstates of the Hamiltonian are states with a fixed to- 
tal spin S and its projection S z ; for a constant S the state 
with the minimum energy correspond to the maximum of 
the projection S z , S z — S. The energy can be written as 
E(S,S Z ) = JS{S+l)/2-gn B HS z , and S z changes dis- 
continuously with jumps from S z = n—l to S z = n at the 
field H n = Jn/(gfJ,B)', The saturated state is reached at 
H = H ex , where the exchange field is H ex = 2JS/ (g^s)- 
Other exactly solvable models for clusters with AFM in- 
teraction are discussed in the review £ In the presence of 
magnetic anisotropy, the quantum problem has no exact 
solution, excluding Ising model that corresponds to ft = 1 
and B = in the Hamiltonian (JTJ. In the latter case, 
which is classical in essence, S z = at H < H ex , and the 
saturated state S z = 2S at H > H ex = JS/(g^B)- 



We study the magnetization curves for a spin pair with 
antiferromagnetic coupling and uniaxial anisotropy in the 
magnetic field directed along the easy axis. It is found 
that the magnetization curve appeared to be more com- 
plicated than for the case of a pure exchange interaction: 
the jumps become nonequidistant and their value AS Z 
can exceed unity. The behavior of the magnetization 
curve also depends on the origin of the anisotropy. It is 
noted similarity of a spin pair and a macroscopic AFM 
close to the point of the spin-flop transition (SFT), that 
it can be useful for the qualitative study of the system. 
We also find out values of the field where discrepancies 
between the quantum and classical model are not small 
even in the limit S — > oo. 



I. QUANTUM MODEL AND ITS 
QUASICLASSICAL ANALYSIS 

Let us consider the Hamiltonian for a pair of spins Si 
and S2 with antiferromagnetic interaction J > 0, uniax- 
ial anisotropy, and an external magnetic field H directed 
along the easy axis z: 



n = j[s z s z 
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Here, both the single-ion anisotropy B and the exchange 
anisotropy nJ are taken into account. We start from 
the classical analysis of the problem. If the operators Si 
and S2 are considered as classical vectors (spin of sub- 
lattices) of a macroscopic sample, a useful resemblance 
with the problem of SFT appears, see Ref ^»i2iii. Denote 
I = (Si — S^VISi — S2I; in the polar parametrization 
l z = cos#, l x + ily = sin#exp(z0). Eliminating the slave 
variable Stat = Si + S2 (the vector of the total spin) 



2g/j,BH(e z — I ■ cos8 
J(2-k + /3cos26>) 



(2) 



where f3 = B / J and e z is the unit vector along the sym- 
metry axis, see details in Reftii, the energy W(9) can be 
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written as 



W(9) = - 



(gfi B H) 2 sin 2 9 
J(2-k + /3cos26») 



+ J{n + (3)S 2 sui 2 9 . (3) 



In the absence of the magnetic field a standard re- 
sult is renewed: both exchange and single-ion anisotropy 
make additive contributions in the effective magnetic 
anisotropy [k + [3)JS 2 sin 2 9 that is easy-axis one at 
k + (3 > 0. Minimizing this energy at H ^ 0, we ob- 
tain various magnetic phases. The collinear phase $|| is 
stable at H < Hi and S z = (9 = or tt) in it. The 
spin- flop phase $x is stable at H > H 2 . ft is character- 
ized by 9 — ±7r/2, and the projection of the total spin 
linearly depends on the field, S z — 2SH/H ex . When 
H > H ex = JS(2 — Hi — /3)/(g/j,B) the saturated state 
with S z = 2S exists. Characteristic fields separated the 
regions of stability are determined as 



JS 



#i = — V(2 -« + /?)(« + 0), H 2 = ^, (4) 

TQ 

H sf = V(2-k-13)(k + P) . (5) 

It is note that in the absence of the single-ion anisotropy 
the expression for the exchange field coincides with the 
saturation field for the exactly solvable quantum models 
of the Heisenberg and Ising magnets. 

If the anisotropy is purely exchange, (3 — 0, the degen- 
erate case with the equal fields Hi and H 2 appears. If 
(3 > 0, the transition from the collinear state $n to the 
spin-flop state <&j_ is discontinuous (SFT is a first-order 
transition). The respective field H = H s f is introduced 
in the Eqs. (0J. The value S z is changed on AS^, where 



AS z cl = 2S* 



(6) 



It is possible another case Hi < H 2 , when the tran- 
sition <I>|| <-> $^ passes via two second-order transitions 
with a skew phase $ / between themiiSiii The model © 
require the competitive anisotropics — the single-ion one 
with (3 < 0, that at k = lead to easy-plane anisotropy, 
and a large exchange anisotropy n > —(3. The effec- 
tive anisotropy can be easy-axis. In the range of fields 
Hi < H < H 2 the value S z changes linearly from to 
AS z \h=h 2 ~ AS^, but the slope of the straight line is 
deeper than in the case of the spin-flop phase, see Fig.H 



II. QUANTUM PROPERTIES 

These results can be juxtaposed with the magnetiza- 
tion curves obtained from the quantum model in the 
limit case of the large spin if wc adopt as a hypothe- 
sis that all linear dependencies S Z (H), see Fig. ^ are 
changed to the step-like ones with equidistant unit jumps 
AS Z = 1, and the jump AS z t > 1 at the point of SFT 
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FIG. 1: The ground state of the classical AFM with single- 
ion and exchange anisotropies for different values of the field 
(schematically). The arrows are spin directions in the differ- 
ent states of AFM. (a) — the single-ion anisotropy prevails; 
(b) — the exchange anisotropy prevails. 



holds in the quantum case. This is reasonable when 
the field of SFT is stronger than the field of the first 
quantum jump H s f > J/(gfiB) for the isotropic model, 
and AS Z [ > 1. Both conditions lead to the inequality 
k + (3 > 1/(2S 2 ). It is valid in the limit case S > 1 even 
for small anisotropy. 

It is convenient to depict the behavior of the system 
in the (b, (3) plane, where b = k + (3 is the constant of 
the effective anisotropy, see Eq. (0. According to the 
classical consideration, SFT with the jump AS z t > 1 can 
occurs only for (3 > 0, and regions with AS z t = n > 
1 are vertical stripes bounded from the bottom by the 
straight line (3 = and from the sides by the vertical 
lines b — b n and b — b n +i, where b n — 8n 2 /(n 2 + AS 2 ). 
The behavior with equidistant unit jumps similar to the 
isotropic model can also be expected for the cases of a 
small anisotropy or an easy-plane anisotropy b < 0. For 
large values of anisotropy b and (3 < the value S z = 
persists up to the field Hi, after that the concentration 
of ri jumps takes place in the narrow range Hi < H < 
H 2 . Thus, the classical approach gives a rich variety of 
types of the magnetization curves. This prediction is 
corroborated by the numerical results obtained for the 
quantum model , but an important exception must be 
made. 
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FIG. 2: Regions with a different behavior of the magnetiza- 
tion curve in the plane (« + /3, 0) for 2S 1 = 10. The lines 
separate the regions of the plane with different values of the 
first jump A5* z (indicated by digits). Regions with even AS Z 
are shaded. The jump with AS Z = 10 occurs when K + /3 > 1; 
the only point (the Ising point k = 1, j3 = 0) correspond it in 
the figure. Letters in circles mark the characteristic regions of 
the phase diagram; the behavior of the magnetization S z (H) 
for them is presented in Fig [3] Bold curve separates region 
where AS Z > 1. The curve below it limits the region with 
lattice of parallelograms where the first jump is AS Z = 1, but 
among subsequent jumps the value AS Z > 1 is present, see 
Fig- C2d) and explanation in the text. 



III. NUMERICAL SIMULATION 

The Hamiltonian commutes with the operator of 
the total spin projection S z = S* + 5|, and the complete 
Hilbert space of the problem decomposes into subspaces 
with fixed values of S z . The operator projected onto 
these subspaces has the form of three-diagonal matrices 
whose eigenvalues can be numerically found using the 
QR-algorithm. Another important simplification of the 
problem is possible because the Zeeman term commutes 
with S z , and the eigenvalues of the Hamiltonian E(S Z , 0) 
obtained for the case H = can be used for any field 
through the shift, E{S Z ,H) = E{S z ,0) - g^ B HS z . 
These properties permit analyze the problem in the com- 
plete phase plane («,/?), even for a sufficiently large val- 
ues 25 ~ 250 in a reasonable time on a personal com- 
puter. 

Numerical investigation shows that almost all proper- 
ties of the classical SFT hold in the quantum model Q), 
even for a small value of spin S ~ 5. First, the state with 
S z = exists in the range < H < H±. Second, jumps 
with AS Z > 1 are also observed for (3 > 0, when the 
classical SFT is a first-order transition. Third, the con- 
centration of jumps in the narrow range Hi < H < Hj 
is observed. The shape of regions with a fixed value 
AS Z > 1 qualitatively corresponds to the classical cal- 
culations even for small spins. Particularly, the curves 
separated regions with different AS Z — n > 1 at (3 > 
are approximatively vertical, see Fig. [2 But qualitative 
discrepancies are also evident. First of all, for small AS Z 
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FIG. 3: Dependence S z on the field for the four characteristic 
regions denoted by letters in Fig. Data is obtained for 
2S — 10. (a) — almost exchange behavior: b = 0.1, (3 = —0.2; 
(b) — merging of jumps, an analogue of SFT: 6 = 0.515, 
(3 — 0.1; (c) — strong concentration of jumps, an analogue 
of a skew phase $z: b = 0.56, f3 = —0.1; (d) — nonunit 
jump AS Z > 1 from state with S z ^ 0, an analogue of the 
phase transition $z «-» $±: f3 = —0.09, b = 0.7, see detailed 
description in Fig.|3 
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FIG. 4: Curves separating the spin-flop and skew phases for 
three values of the spin: o - 25* = 8; A - 2S = 16; • - 2S = 32; 
the continuous line denotes data for 25* = 256. The dashed 
rectangle is the region where the minimum of the continuous 
line occurs. In the inset — this region and the part of the 
curve adjacent to the j/-axis for 2S = 256. 



regions with A5 Z > 1 extend into the lower half-plane. 

The behavior of the line separating regions with AS Z = 
1 and A5 Z > 1 is astonished. The discrepancy from the 
classical result (3 = is not small even for large spins 
25 ~ 200, see Fig. El The extent of this region in the 
direction of (3 grows when the spin increases and arrives 
at the limit value (3 ~ —0.3, but its area decreases with 
5 — > oo. 

For magnetization curves of the kind d, see Figs. Eland 
EJd), in the skew phase nonunit jumps are observed, that 
corresponds to the first-order transition $^ «-> $x m the 
limit case of large spins, but the classical theory predicts 
that this transition has a second order, and the region 
d is absent. Also the nonequidistantness of jumps with 
AS* 2 = 1 in the region of the concentration of jumps (be- 
low the bold line in Fig. El is observed, see Fig. EJc). 
Their merging is illustrated in Fig. Efd). I n the latter 
case that corresponds to the narrow region d in Fig. El 
the magnetization process is the following: several first 
jumps have AS Z = 1, then the jump with AS Z > 1 pro- 
ceeds (always one), and after that equidistant jumps with 
AS Z = 1 appear. In terms of the spin-flop analogue this 
behavior can be treated as that the transition $^ <-> $^ 
is of the first order, which is impossible for the model 
given by the expression ©■ The more detailed struc- 
ture of the region d is presented in Fig. El for 25 = 64. 
Parallelograms in the phase diagram (6, (3) , the region d 
in Fig. El corresponds to transitions from a state with 
S z ^ to state with S z + AS Z , AS Z > 1. Parallelo- 
grams with S z and S z ± 1 is adjoint, and parallelograms 
with equals AS Z are attached at corners. This structure 
is responsible for the sawtooth shape of curves separat- 
ing regions with AS Z > const in Fig. El However, the 
isolines in this figure and the plot in the inset show that 
the value of the jump decreases rapidly removing from 
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FIG. 5: Curves separating regions with a fixed jump value 
in the transition <E>z «-» <&x for 2S — 64. Above the bold 
line the jump AS Z > 1 occurs from a state with S z = 0, see 
correspondent lines in Figs.|5]and|l] Below this line there are 
regions where the nonunit jump AS 12 denoted by inequality 
occurs from a state with S z > 0. The section of the main 
plot for k + P — 0.5 (dotted line) is shown in the inset. The 
axis y is the logarithm of the value of the jump normed on 
ASS = AS Z | M . 



the line of separation states with AS Z 
in the direction of negative (3 values. 



1 and AS Z > 1 



IV. RESULTS AND DISCUSSION 

Thus, we find out a complicated behavior of the step- 
like magnetization curves for a pair of quantum spins 
with AFM interaction and easy-plane anisotropy that de- 
pends on the two anisotropy constants with different ori- 
gin. If the single- ion anisotropy prevails, "merging" of 
unit jumps into the jump with AS Z > 1 can take place. 
In the opposite case, when the exchange anisotropy dom- 
inates, all jumps are unit, but they are concentrated in 
a narrow field range. We reveal a qualitative similarity 
of the quantum problem and the corresponding classical 
one down to S ~ 5. These observations can be suitable as 
a starting point for the analysis of the experimental data 
ans studying the influence of other kinds of interactions, 
but it is far beyond this brief paper. At the same time, 
there are quantitative discrepancies between the classi- 
cal system and the quantum one. These discrepancies 
are important in the range of the fields where the transi- 
tion from AS Z = 1 to 2 occurs. The quantum effects that 
are free of a small value 1/5 manifest themselves there. 
Apparently, the explanation of this unexpected fact is 
associated with the follows: for the small anisotropy the 
quantum model has a singlet ground state that is strongly 
dissemble from the quasiclassical Neel picture. Thus, the 
region of a pure quantum behavior is revealed for the 
system in the limit 5 — > oo. Such effects are possible in 
the physics of antiferromagnetism. We mention the Hal- 
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dane hypothesis 12 that the ground state of a quantum spin pair considered above demonstrates that antiferro- 
spin chain is different for integer and half-integer spin magnets can manifests substantially quantum properties 
even in the limit S ^ oo. Thus, the simple model of a even for large spins. 
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